Velocity profiles in shear-banding wormlike micelles 
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Using Dynamic Light Scattering in heterodyne mode, we measure velocity profiles in a much 
studied system of wormlike micelles (CPCl/NaSal) known to exhibit both shear-banding and stress 
plateau behavior. Our data provide evidence for the simplest shear-banding scenario, according to 
which the effective viscosity drop in the system is due to the nucleation and growth of a highly 
sheared band in the gap, whose thickness linearly increases with the imposed shear rate. We discuss 
various details of the velocity profiles in all the regions of the flow curve and emphasize on the 
complex, non-Newtonian nature of the flow in the highly sheared band. 

PACS numbers: 83.60.-a, 83.80.Qr, 47.50.+d, 83.85.Ei 



Understanding the correlation between mechanical and 
structural response in non-Newtonian fluids submitted 
to high deformation rates is crucial on both fondamcn- 
tal and technological grounds jj]]. Among the variety 
of complex fluids investigated in recent years, a wide 
class exhibits flow-structure coupling that leads to a 
strong shear-thinning behavior: along the steady-state 
flow curve (shear stress a vs. shear rate 7), a drop of 
up to three orders of magnitude in the effective viscos- 
ity r\ = a/j is observed in a very narrow stress range 
leading to a stress plateau (for a review, see for instance 
Refs. jl], In correlation with this stress plateau, bands 
of different micro-structures and normal to the velocity 
gradient appear. Such bands correspond to a new shear- 
induced structure (SIS), whose low viscosity is in gen- 
eral supposed to be responsible for the shear-thinning. 
This so-called shear-banding behavior has been observed 
in both ordered mesophases (lamellar, hexagonal, cubic) 
H and transient gels Q]. 

A particularly well-documented example is the group 
of wormlike miccllar systems of self-assembled surfactant 
molecules H, ^ . They consist of very long cylindrical ag- 
gregates whose configurations mimic polymer solutions. 
However their dynamics is strongly modified by the equi- 
librium character of the chains, which enables them to 
break and recombine fij] . Generically, one starts from an 
isotropic viscoelastic solution of these micelles above the 
scmidilutc regime, which behaves like a Maxwell fluid at 
low shear rates. Upon increasing 7 and entering the non- 
linear regime, the onset of the stress plateau for a critical 
shear rate j t is associated with the nucleation and growth 
of highly bircfringent bands, suggesting strong alignment 
of the micelles along the velocity direction @, . As the 
shear rate is further increased above j 1 , the new organiza- 
tion progressively fills the gap at almost constant stress, 
up to a second critical shear rate % . Above 72, the sys- 
tem enters a second regime of apparently homogeneous 



structure, with a second branch of increasing stress. The 
flow curve of Fig. [j] is typical of miccllar systems like that 
investigated in the present work. Such a stress plateau 
has been reported for concentrations close to the equilib- 
rium isotropic ncmatic (I-N) transition, where coupling 
between the order parameter and the flow could explain 
an out-of-equilibrium phase transition to local nematic 
order || |8|. This behavior has also been observed in 
more dilute systems, where arguments based on insta- 
bilities in the underlying flow curve have been invoked 
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FIG. 1: Steady-state flow curve for a 6 % wt. CPCl/NaSal 
solution in 0.5 M brine at 21.5°C. The data were collected 
under controlled shear rate in a 1 mm gap Couette cell. Ver- 
tical lines indicate the limits of the banding domain. For 
7 > 72 ~ 26 s _1 , the data were fitted by a power-law 
fluid a = 35.9 7 019 (dashed line) and by a Bingham fluid 
— 55.1 + 0.43-y (solid line). Inset: high shear branch; error 
bars account for temporal fluctuations of a. 

In the earliest simple pictures of the flow along the 
stress plateau, the system is supposed to separate into 
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two differently sheared bands: a weakly sheared region 
that flows at j ± and a highly sheared region at 7 2 ||, |). 
Experimentally, this very basic issue of the nature of the 
flow field remains rather obscure. Among previous stud- 
ies, Nuclear Magnetic Resonance (NMR) imaging has 
proved the most successful technique to measure local ve- 
locity and/or local shear rate. Callaghan and coworkers 
have shown for the first time the existence of inhomo- 
geneous velocity fields in wormlike micelles O, [l2| . 
However, due to the limited amount of data in NRM 
images, no clear evidence supporting the simple shear- 
banding scenario has been provided. In particular, one 
of the most detailed NMR studies promotes a quite dif- 
ferent dynamical description: instead of being submitted 
to a high shear rate, the nucleated nematic-like structure 
seems to undergo a zero-shear-rate plug-like flow and to 
move like a piece of gel, probably involving important 
wall slip and fracture behaviors jl2). This last observa- 
tion has put the whole naive shear-thinning picture on 
uncertain ground. 

In this Letter, we show that for the much studied 
wormlike micellar system CPCl/NaSal in brine and for 
concentrations far from the I-N transition, the simplest 
scenario holds. The originality of our work relies on 
recording both the local velocity and global rheological 
data in Couette geometry simultaneously along the whole 
flow curve. This enables to demonstrate for the first time 
the nucleation of a highly sheared band at a critical stress 
and to follow its growth from the rotor to the stator as 
the shear rate spans the stress plateau. We show that 
the width of the band grows linearly with the shear rate, 
in agreement with the decrease of the measured effective 
viscosity, and that the SIS is not Newtonian. 

Our system of elongated wormlike micelles consists of 
a binary mixture of cetylpyridinium chloride (CP + , Cl~) 
and sodium salicylate (Na + , Sal - ) in 0.5 M NaCl-brine. 
Since the pioneering work by Hoffmann and Rehage [jl3| , 
this extensively studied system has been demonstrated to 
exhibit a stress plateau and optical birefringence shear- 
banding We focus on a 6 % wt. sample at 21.5°C, in 
the domain above the scmidilutc regime (0.5-5 % wt.) 
but far below the equilibrium I-N transition (around 
20 % wt.). For the concentrations considered here, the 
system is known to be a perfect Maxwell fluid in the lin- 
ear regime, with a typical relaxation time of about 1 s. 

Rheological flow curves are measured using a standard 
rheometcr (TA Instruments AR 1000) and transparent 
Couette cells of outer radius 25 mm and different gap 
widths e (e = 1 or 3 mm). The temperature is main- 
tained at 21.5±0.1°C by a water circulation around the 
cell. To access local velocity, we use a heterodyne Dy- 
namic Light Scattering (DLS) technique that has been 
described elsewhere |l4|]. The measurement of the local 
velocity relies on performing the interference between a 
reference beam and light scattered from a small volume 
of the sample of typical size 50 /im. The correlation func- 



tion of the interference signal exhibits oscillations at the 
Dopplcr frequency qv where q is the scattering wavevec- 
tor and v the local velocity. Good statistical convergence 
is achieved by averaging the correlation function over 3 s. 
Velocity profiles are measured by moving the rheometcr 
along the velocity gradient by steps of 30 /im. This tech- 
nique enables us to obtain a complete velocity profile in 
about 2 min, more than 10 times faster than NMR ve- 
locity imaging. Typical uncertainties are about 5 %. 

In order to enhance the scattering properties of our 
system, we add a small amount (1-5 % wt.) of 30 nm col- 
loidal particles (Ludox from Aldrich). We checked that 
the rheological properties, particularly the plateau be- 
havior, was not affected by the addition of those scat- 
terers. Indeed, the small size of the scattcrers compared 
to the typical mesh size of the micellar network in this 
concentration domain should lead to a negligible influ- 
ence on the structural and mechanical behaviors of the 
sample. 
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FIG. 2: Velocity profiles recorded while measuring the flow 
curve of Fig. |l[ (a) j = 1 (o),5 (•), and 12 s" 1 (■). (b) 
7 = 15 (o),22 (•), and 28 s' 1 (■). The value x = (x = 1 
resp.) corresponds to the rotor (stator resp.). The dotted 
line is v(x) = ji (e — x) with 71 = 2.5 s _1 . The dashed lines 
correspond to a shear rate 72 = 26 s -1 



Figure H shows the steady-state shear stress recorded in 
our micellar system at various imposed shear rates. This 
flow curve presents a stress plateau at a « 65 Pa that 
extends from j ± « 2.5 s _1 to % ~ 26 s _1 corresponding 
to a drop in the effective viscosity by a factor of 10. Note 
that for 7 > 7 2 , the stress response is no longer strictly 
stationary: a fluctuates by about 2 % around its mean 
value. Finally, at 7 « 37 s _1 , the sample began to frac- 
ture so that no measurement is available at higher shear 
rates. 

Figure ^ summarizes our main result: velocity pro- 
files in the gap are displayed for 7 ranging from 1 s _1 
to 28 s _1 . For all the shear rates in the plateau do- 
main, the velocity profiles exhibit an unambiguous band- 
ing structure: two regions of different well-defined local 
shear rates coexist when j 1 < 7 < 72 . The highly sheared 
band nucleates on the inner wall of the Couette cell at 
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the onset of the plateau and expands as 7 is increased 
up to 7 2 where the band fills the whole gap. As can be 
seen on Fig. 0, the local shear rate in the weakly sheared 
region remains constant and equal to A f 1 = 2.5 ± 0.2 s" 1 . 
Moreover, although significantly curved (see below), the 
velocity profiles in the highly sheared region are compat- 
ible with a local shear rate of 72 = 26 ± 1 s _1 . Note that 
in all our data, no noticeable wall slip was observed. Fi- 
nally, these profiles easily yield the width h of the highly 
sheared band. Figure || demonstrates that h increases 
linearly with 7. This corresponds to the simplest rela- 
tion between the imposed shear rate and the two local 
values j 1 and j 2 : 7 = (1 — a)j 1 + 073, where a = h/e 
increases from to 1 over the plateau region. 
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FIG. 3: Width h of the highly sheared band as a function 
of 7. With the value 7i = 2.5 ± 0.2 s _1 inferred from the 
velocity profiles of Fig. ti, a good linear approximation of h(j) 
is obtained with 72 = 26 ± 1 s _1 (solid line). 



FIG. 4: Velocity profiles measured in the 3 mm gap Couette 
cell and rescaled by the rotor velocity v for two values of 7 
below the banding transition. For 7 = 1 s _1 (o), the nor- 
malized data are well fitted by a weakly shear-thinning law 
a ^, 7 ±°- 1 (dotted line). For 7 = 2 s _1 (•), a stronger 
shear-thinning behavior a ~ ^,o.28±o,03 pj-Qyjdes a very good 
fit (dashed line) . The solid line is the velocity profile expected 
for a Newtonian fluid. 
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Beyond shear-banding evidence, the resolution of our 
technique enables us to analyze in more detail the flow 
behavior of the material in the regions below and above 
the plateau. In Fig. |J are displayed two normalized veloc- 
ity profiles below the plateau at 7 = 1 s _1 and 7 = 2 s -1 
in a 3 mm gap. For 7 = 1 s , the velocity profile is 
very close to a straight line, consistent with the Newto- 
nian behavior of the miccllar solution at low shear rates. 
However, as shown in the inset of Fig. [|, our data do not 
exactly fall on the Newtonian velocity profile but rather 
present a small curvature which can be accounted for by 
a power-law a ~ ,y0.7±o.i_ when 7 = 2 s _1 i.e. just 
below the onset of shear-banding at 7 X , the curvature is 
much more pronounced and the power law a ~ <y°- 28 ± - 03 
yields a perfect fit of the data. This demonstrates the ex- 
istence of weak shear-thinning on the low shear branch, 
that sharply increases as the banding transition is ap- 
proached. 

Looking more closely at the velocity profiles in the 
plateau region, one notes that the highly sheared band 
displays a significant curvature even in a 1 mm gap (see 



FIG. 5: Velocity profiles in the high shear branch at (a) 7 = 
28 s _1 and a = 66.9 Pa and (b) 7 = 32 s _1 and a = 68.5 Pa. 
In each case, two solid lines correspond to two different but 
undistinguishable fits: (i) a power-law behavior a — 35.9 7 019 
and (ii) a Bingham behavior a = 55.1 + 0.43 7. The dotted 
lines account for a global stress fluctuation of 2 %. 

Fig. ||(b)). This shows that the SIS cannot be described 
by a Newtonian fluid. The same observation holds for 
velocity profiles recorded at 7 > 72 where the SIS has 
completely invaded the gap. A precise analysis of those 
velocity profiles may help to understand both the global 
rheological measurements and the behavior of the SIS. 
Indeed, if one attempts to fit the flow curve of Fig. for 
7 > 7, according to a power-law a = Aj n i the result 
is dramatically shear-thinning: a = 35.9 7 019 (see inset 
of Fig. |l|). Such a shear-thinning exponent n allows us 
to account for the curvature of all the velocity profiles 
at 7 > 7 2 with a single prefactor A consistent with the 
rheological data (sec Fig. || for the cases 7 = 28 and 
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32 s^ 1 where the stress has been allowed to vary by 2 % 
according to the measured temporal fluctuations of a) . 

However, above %, other constitutive relations cr(j) 
may fit the experimental data as well. For instance, a 
Bingham fluid a = a + Aj with a = 55.1 Pa and 
A = 0.43 Pa.s also closely fits both the flow curve and ve- 
locity profiles for 7 > 7 2 with a single set of parameters 
{a ,A) (see inset of Fig. [l] and Fig. ||). Such a Bing- 
ham behavior is suggested from the analogy with other 
related viscoelastic systems, whose similar nonlinear be- 
havior has been interpreted in terms of a steady-state 
system of bulk fractures Here, the range of 7 on the 
high shear branch is limited to 26-37 s _1 so that we can- 
not discriminate between the two above behaviors and 
only further studies will help to select the correct non- 
Newtonian relation for the high shear branch. 

Let us now further discuss our results in light of pre- 
vious studies. Among the great diversity of wormlikc 
micellar systems, a universal feature seems to emerge: 
the presence of a robust stress plateau associated to the 
growth of a SIS. However the present results and recent 
local NMR velocity measurements by Fischer et al. JlJ] 
lead to two radically different descriptions of the flow 
field. Ref. |l2j suggests the nuclcation and growth of 
a gel whereas our data clearly shows the coexistence of 
two different well-defined shear bands flowing at jt and 
7 2 respectively. 

These important discrepancies may be explained by 
two different arguments. First, the distance to the I- 
N transition is drastically different in the systems under 
study. Indeed, Fischer et al. used a CTAB/D 2 system 
at a surfactant concentration of 20 % wt. very close to 
the I-N transition (w 21 % wt. at the working tempera- 
ture). In our case (CPCl/NaSal in brine), the surfactant 
concentration of 6 % wt. was much more distant from 
the I-N transition. Second, significant temporal fluctua- 
tions of the flow field are mentioned in both the present 
work and Refs. p). If the flow of the SIS is time- 
dependent, we think that the data analysis in the two 
techniques may be affected in different ways. The tem- 
poral resolution of our heterodyne DLS setup provides an 
accurate value of the local velocity averaged over a few 
seconds. With NMR imaging, a broad velocity distribu- 
tion is reported in the SIS, sometimes showing several 
peaks. We believe that this rich spectral content may 
complicate the interpretation of the average NMR data. 

Another important open issue concerns the very na- 
ture of the SIS and its rheological behavior. As far as 
the structure is concerned, birefringence appears as a ro- 
bust experimental fact || nj. However, the correlation 
between birefringence and local nematic order is only 
confirmed for systems close to the I-N transition JlT| . 
This points to possible differences between shear-induced 
structures far from and close to this equilibrium transi- 



tion. As for the behavior under shear, flow instabilities in 
the SIS at high shear rates leads to fracture and ejection 
of the sample from the rheometer || |j| . We believe that 
the non-Newtonian features of the SIS revealed in this 
Letter may play an important role in such instabilities. 
Indeed temporal fluctuations of the stress and of the bire- 
fringent bands seem to suggest intermittent and localized 
fracture-like flow events, yet not directly observed so far. 

Although separate studies of rhcology, birefringence, 
and more recently NMR velocity profiles, have been per- 
formed, a complete picture of the nonlinear rhcology of 
wormlike micelles remained rather elusive in the litera- 
ture due to the lack of experimental data. Coupled with 
optical birefringence and rheological measurements, our 
heterodyne DLS setup appears as a well-suited tool to 
provide the missing information on this much studied 
problem. 
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